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ABSTRACT 

An investigation into the limitations on the enhanced field-emitted current density 
in a fast-pulsed (rise-time = ns), high voltage (> 10 6 V), 1-inch vacuum diode was 
conducted using a computer simulation based on the Fowler-Nordheim equation. 

Oscillations in the emitted current density (due to the change in the amount of 
space charge within the gap) were found to quickly decay into a final steady-state for 
the voltages applied. Steady-state values for a wide variety of work functions, electric 
field enhancement factors (based on the theory that "whiskers" on the cathode surface 
experience varying degrees of enhancement), and applied potentials were compared 
to two benchmarks: the amount of current density required to explode a whisker in 
<10 ns by joule heating ( J E = 10 9 A/cm 2 ); and the Child-Langmuir (C-L) space- 
charge-limited current density. 

Steady-state values were found to be less than J E . One model of the formation 
process of a plasma at the cathode surface requires that J E be met or exceeded by 
the steady-state value. Thus, such a model is not supported by this project’s findings. 

The C-L limit is based on a thermionic-type emission process. As only pure field 
emission (i.e., no thermionic emission included) was considered, the steady-state 
values were, in all conclusive cases, less than the corresponding C-L limited values. 
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I. INTRODUCTION 



The breakdown of the surface of the cathode in a pulsed vacuum diode has been 
the subject of much investigation in recent years. The existence of crater-like cathode 
spots after breakdown can be explained by unipolar arcing, but the phenomena 
occurring prior to the arcing are still debated. 

One model contends that a metallic "whisker" on the surface of the cathode 
provides a source for the enhanced emission of electrons. After a period of time (less 
than 10 ns for a fast-pulsed diode), the whisker explodes (due to the rapid build up of 
energy produced by the emission process) into a dense plasma which in turn drives 
the unipolar arc. 

This project will attempt to show that field emission alone (i.e., not including 
thermionic emission) from a whisker cannot reach high enough values to trigger the 
explosion mechanism. A computer simulation utilizing the Fowler-Nordheim equation 
for field emission will be used. The results will also be compared to the well-known 
Child-Langmuir space-charge-limited current density for a vacuum diode. 

A look at the dimensionless variables used in the program and the program itself 
are included in the appendices. 
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II. BREAKDOWN AND EMISSION THEORY 



A. BACKGROUND 

For over six decades scientists have been investigating pulsed-voltage vacuum 
diode breakdown. With the advent of newer technologies, shorter and shorter 
voltage-pulse rise times have become possible, with values ranging from 10 -6 seconds 
in 1928 to 10 -9 seconds today [Ref. 1]. Corresponding to today’s very short rise times 
has been the discovery that breakdown even occurs on the nanosecond time scale. 
The effects of the breakdown process on the surface of the cathode have been well 
documented [Refs. 1-3], and a theory of the formation of cathode spots (Fig. 1) via 
unipolar arcing 1 has recently been proposed [Refs. 2,3]. In dispute, however, is the 
process by which the plasma driving the unipolar arc is formed. 

Cathode spots are characterized by the formation of a dense quasi-neutral 
plasma near the cathode surface. There are two basic schools of thought concerning 
the primary mechanism for breakdown and the source of the plasma, only one of 
which will be attended to here. Anode-initiated breakdown would occur when ions 
released from the anode (due to bombardment of the anode by electrons emitted from 
the cathode) interact with the cathode to form the plasma. Here fast-pulse 
( = nanoseconds) breakdown shall be investigated, in which the anode-emitted ions 
would not have time to reach the cathode (for the 1 in = 2.54 cm diode gap of interest) 
before breakdown occurs; anode-initiated breakdown shall be discussed no further. 2 

1 The concept of an arc which originates and terminates on the same surface (i.e., "unipolar") 
was first expounded by Robson and Thonemann in 1958, and was later expanded by Schwirzke 
in the 1980’s [Refs. 2,4], 

2 For a more thorough treatment of anode-initiated breakdown, see [Ref. 1], 



2 



j* r *|*''v-*»**' 



1 

100 ym 






m 

• ♦ 



• « * 9 









1 



* 



#% 



#^(i^ 

V»W* 



Figure. 1. Crater-like cathode spots formed by unipolar arcing. 
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In the cathode-initiated breakdown school, it is proposed that all of the ions and 
electrons which make up the plasma are produced by processes which are initiated by 
events occurring only at the cathode [Ref. 1]. Within this school, however, there are 
many camps, each with its own views on what processes do indeed occur, and which 
events really take place to initiate them. One of the more reknown camps is led by 
Schwirzke, who recently proposed that neutrals released from the surface of the 
cathode are ionized by the emitted electrons, thus forming a plasma near the cathode 
surface [Refs. 2,3]. Another is led by Dyke (et al) 3 who proposed that joule heating of 
an electron emission site causes the site to vaporize, and in the process ions are 
created which allow the plasma to form [Ref. 1]. Here shall be examined a third camp, 
championed by Mesyats [Ref. 1], in which Dyke’s theory is taken to its most extreme 
limit, that of an explosion of the emission site into a plasma directly. The goal is to 
show that an elementary version of a Mesyats-type "explosive" model cannot support 
breakdown on a time-frame indicative of fast-pulsed voltage vacuum diode breakdown. 

B. THE MODEL DIODE 

The model's anode and cathode will consist of two polished metal plates 
separated by a vacuum. Despite the polishing and the degassing effect of the 
vacuum, however, some small-scale "irregularities" will nevertheless exist on the 
surfaces of the plates. Weakly bound bits of dust, oxides, and non-metallic inclusions 
litter the cathode surface, providing excellent sites for the enhanced field emission of 
electrons [Ref. 2]. These so-called "whiskers" will vary in size and shape, with the 
more needle-like (height-to-base ratio » 1) metallic protrusions providing the most 



3 See [Ref. 1] for a list of references on Dyke’s experiments. 
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enhancement of the field emission. Flatter, non-metallic (e.g., dielectrics, semi- 
conductors, etc.) sites can also provide significant enhancement [Ref. 5]. 

C. ELECTRON EMISSION FROM THE CATHODE 

There are two mechanisms for the electric-field-enhanced emission of electrons 
from a metallic surface: thermionic (high temperature, low applied field); and field (low 
temperature, high applied field) [Ref. 4]. The presence of the electric field at the 
surface lowers the potential barrier at the plate-vacuum interface, thus allowing 
electrons in the conduction band an opportunity to become free from the metal (Fig. 2). 
In thermionic emission, electrons at high temperatures gain enough kinetic energy to 
pass over the barrier and into the gap. In field emission, the electrons tunnel through 
the barrier instead. In general, the electric field at the emitter surface must be on the 
order of 10 7 V/cm for field emission to take place [Ref. 6]. For a diode gap of 2.54 cm 
(= 1 in), an applied potential of 2.54x1 0 7 V would be necessary to achieve field 
emission, but as noted previously, the whisker emission sites experience an enhanced 
applied electric field (Fig. 3). The enhancement factor, m, varies from 1 to several 
hundred, with needle-like whiskers having the largest values [Ref. 7). For electric 
fields less than 10 7 V/cm, thermionic emission will dominate for temperatures above 
approximately 1500 K, with little substantial emission at lower temperatures [Ref. 6). 

For the case of the non-metallic inclusions, a combination of field and thermionic 
emission occurs. Electrons tunnel from the metallic cathode material into the 
dielectric, then are thermionically emitted into the gap (see Fig. 4) [Ref. 5]. 
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Figure 3. Enhancement of electric field at whisker tip [Ref. 7]. 
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Figure 4. Combination of field and thermionic emission for a dielectric [Ref. 5]. 



For this elementary model, only field emission from metallic whiskers shall be 
considered, with the following assumptions: 



1. The time-scale for the explosion of a whisker is only a few 

(< 10) nanoseconds; 

2. The applied voltage rises instantaneously to a constant value at t = 0; 

3. The gap is perfectly evacuated; no positive ions are present, so no heating of 

the site by ion bombardment is possible; 

4. The initial temperature of the cathode (=300 K) is too low for any 

appreciable amount of thermionic emission to occur at first; and 

5. Any increase in temperature during the short time-frame prior to explosion 
shall be neglected (and thus so shall thermionic emission) in this elementary 
model. 
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Pure field emission is governed by the Fowler-Nordheim equation [Ref. 6] 



J(F) = ( 1 .54x10 ) V e 

<t> 



a 3/2 

p2 ( -6.83x1 0 7 ) f(y) -2— 
~o \ 1 rs r 



( 1 ) 



Since f(y) (the Nordheim elliptic function) is slowly varying, it can be approximated as 
a constant 4 



f(y) = 0.94 , 



which leads to (see Figs. 5a, b) 



* 3/2 



p2 ( -6.42x1 0 7 ) 

J(F) = ( 1.54X10 -6 ) — e F 



(2) 



where 



J s emitted current density (A/ cm 2 ); 

F = -mE s enhanced electric field 5 at the whisker tip (V/cm ); 
m s electric field enhancement factor; and 
<t> s work function of the cathode metal (eV). 



J is a function only of F for any given whisker; any negative (i.e., electron) space 
charge which exists in the gap will reduce the emitted current density. If the applied 
potential remains constant, it can be predicted that J will start at some initial value, J 0 , 
and decrease as more and more space charge enters the gap. The electrons, once in 
the gap, will be accelerated across to the anode, eventually disappearing there. When 
the first emitted electrons reach the anode, the electric field at the cathode will 
increase due to the loss of space charge. The increased field will cause an increase 

4 For a more thorough treatment of f(y), see [Refs. 4,6], 

5 Eq. (2) requires that F be defined as a positive quantity, hence the negative sign. 
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Figure 5a. Log plot of J vs. F (Eq. 2). 



14 




Figure 5b. Expanded log plot of J vs. F for various work functions (Eq. 2). 
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in J, but soon the increased addition of space charge will exceed the loss at the 
anode; the amount of space charge in the gap will then rise, the electric field at the 
cathode will decrease, and J will fall. This oscillatory action should diminish as the 
emitted current density at the cathode comes to equal that going into the anode. In 
other words, the current density becomes constant at some final space-charge- 
determined value, J F . 

It should be noted that in reality, the applied potential is not constant from t = 0. 
Rather, it ramps up from zero to some final value. The ramping has the effect of 
causing the current density to rise in a similar fashion; oscillations such as those 
described above are not observed. 



D. NON-APPLICABILITY OF CHILD-LANGMUIR LAW 



It is somewhat important to note, at this point, that J F is not the well-known 
Child-Langmuir space-charge-limited current density (A/cm 2 ) 



Jcl 





v 3/2 

"d 2- 



( 3 ) 



with 



£o s permitivity of free space = 8.85x1 0 -12 C 2 /N-m 2 ; 
e ■ I electron charge I = 1.60x10~ 19 C ; 

M h electron mass = 9.1 1x1 CT 31 kg ; 

V s applied voltage (V) ; and 
d s diode gap length (cm) . 

Jp should be, in fact, somewhat less than J CL . This difference arises due to the fact 
that J CL is not derived for the case of a pure field emission mechanism. Rather, there 
must be an ample supply of electrons (eg., in the case of electrons emitted from a 
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heated filament) always available for emission. In the derivation of J CL , one of the 
boundary conditions is that the electric field at the cathode go to zero due to the 
presence of the space charge (see [F.efs. 4,6,7]). If this were to happen to a pure field 
emitter, then J would go to zero, and a steady state would not develop. Since the 
field emitter cannot emit at J CL for any appreciable length of time without shutting itself 
off, its steady state emission value must be something less: 

J F < Jcl ( P ure field emission ) . (4) 

E. EXPLOSIVE MODEL CURRENT DENSITY REQUIREMENT 

It has been shown by Schwirzke [Ref. 2] that, for stainless steel, a current density 
of J E = 10 9 A/cm 2 would be required to explode a whisker on the cathode surface in a 
few (<10) ns by joule heating. Thus, a whisker would need to emit a J av g = Je f° r 
approximately 10 ns before exploding. If the transient oscillation period is small 
compared to 10 ns, then J avg = J F , and 

J F > J E (explosive emission) (5) 

is required for the explosive emission model. This project will attempt to show that, for 
a wide range of applied potentials, work functions, and enhancement factors, Eq. (5) 
cannot be satisfied, thus proving that field emission alone is not capable of providing 
the current density required by the explosive model. 
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III. SIMULATION 



A. BACKGROUND 

The type of simulation model used here is, as [Ref. 8] so aptly puts it, one of the 
"simplest" possible, viable since the first digital computers came into use in the 1950’s. 
While almost primitive in comparison to the latest innovations in so-called "particle-in- 
cell" simulations [Ref. 8], the model used here nevertheless is time-tested and has 
worked quite well for the rather simple events under investigation. Although the basic 
techniques are over three and a half decades old, the actual program is original and is 
included (in one form) in Appendix B for scrutiny. The program is written in the "C" 
language and was run on a Sun SPARCstation 2 incorporated into a system running 
SunOS™ 4.1.1.. 

B. SIMULATION THEORY 

In general terms, the simulation works as follows. Space charge is simulated not 
as a continuous beam of electrons, but rather as a series of infinitely thin, numerous 
charged disks of varying charge density o [C/m 2 ]. In each time interval, the forces on 
the disks within the gap are computed, the disks are moved, and a new disk is created 
at the cathode. When a disk reaches the anode, it is no longer used in any further 
calculations. In the final steady state, all of the disks will have the same charge 
density; for every time step just as much charge is lost at the anode as is created at 
the cathode. 
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1. Electric Field 



The charge density o is calculated by multiplying the time step At by the 
current density J (Eq. 2) 6 

a = -J At . (6) 

The field-emitted current density at the cathode surface depends on E^O), 
the value of the electric field at the cathode due to all of the space charge disks, and 
Egpp , the component of the electric field at the cathode surface due to the applied 
potential across the diode. In order to find E sc (0), one starts with the one-dimensional 
Poisson's Equation [Ref. 6] 

d 2 V(x) _ q( n i ~ n e ) 

dx 2 eo 1 

which, for the ion density n, = 0 and the electron charge q = -e becomes 

d 2 m = en ( x ) t (8 ) 

dx 2 Eo 

where 

V = potential (V) ; 

x = position ( = 0 at cathode, = d at anode ) (m) ; 
e = I electron charge I = 1.60x10" 19 C ; 
n = electron number density (m -3 ) ; and 
eo = permitivity of free space = 8.85x1 0" 1 2 C 2 /N-m 2 . 



6 The current density J will always be defined as positive, hence a negative sign has been ad- 

ded on the right side of Eq. (6). 
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Using dV(x)/dx = -E(x), one gets 



dE(x) -e n(x) (Q) 

dx £o 

Integrating both sides over the gap, using the relation E sc (d) = -E sc (0), and 
rearranging yields 

E s C (°) = 2^ J 0 e n M dx • 0°) 



Now, using the charged disk concept, the integral becomes a summation 



Esc(O) — n Z e n k Ax . 

^^0 gap 

where n k s electron number density [m -3 ] in disk number k. 
e n A x = -o yields 



( 11 ) 

Realizing that 



Finally, using Eq. (6), 



Esc(O) 



1 

2eo 



Z ~°k • 
gap 



E S c (0) 



1 

2eo 



Z JkAt . 

gap 



(12) 



(13) 



For a disk located at a point x' in the gap, E sc (x') can be determined from 
Eqs. (9) and (10) by breaking the integral of the left-hand-side into two parts and 
subtracting: 

j dEsc(x) = 7 -J 0 -e n(x) dx (14a) 

Esc(O) ®0 J ° 

becomes 

E sc (x')-E sc (0) = X J k At , (14b) 

EO 0-*x' 
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and 



-sc ^ 

J dE sc (x) = J- J _en(x)dx (15a) 

E.cV) * * 

becomes 

Esc(d) - EscOO = - 7 - £ Jk At . (15b) 

^ x' — *d 

Subtracting Eq. (15b) from Eq. (14b) and rearranging yields 

E S c OO - -x- 2 J k AI - -L I J k At . (16) 

^£0 x'->d «0 0->x' 

Eq. (16) is equivalent to saying that the space charge "in front" (x'-»d) will contribute 
positively to the electric field at x', whereas the space charge "behind" (0->x') will 
contribute negatively. 

If the applied voltage remains constant, so will the applied electric field E app . 
The total electric field at x' is therefore 7 

Ej(x') = -E app + Z 4 At - I JkAt . (17) 

^£o x'-+d ^£o o~*x' 

At the cathode, the enhancement factor must be included: 



F = -m E t ( 0) = m J E 



app 2 ea ^ A 1 ' 
gap 



(18) 



A schematic diagram of the electric field components in the diode is shown in Fig. 6. 



7 Like F, E w will be defined as a positive quantity, so the appropriate sign change must be 
made. 
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Figure 6. Schematic diagram of electric field components in diode. 
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2. Equation Of Motion 

The nature of the large (> 10 6 V) applied potentials used in the simulation 
requires that relativistic equations of motion be used. For any electron in a disk at a 
position x, the relation is (based on mass x acceleration = charge x E-field) 

^ ■ jk ■ Tk ■ (19) 



where 

Y = 1 / Vi - v 2 / c 2 ; 
j? s V/ c ; 

t = time (s) ; 

e s I electron charge I = 1 .6x1 0 -19 C ; 

M s electron mass = 9-llxlO -31 kg ; 
c s speed of light in vacuum = 3.0x10 s m/s ; 

(x) = total electric field at position x (V/m) ; 
x = position ( = 0 at cathode, = d at anode ) (m) ; and 

V s velocity (m/s) . 



Substituting into Eq. (19) and carrying out the derivative gives 




Taking the derivative of y and substituting in yields 

With a little more substituting, Eq. (21) can be solved for d\f / dt: 



dV 



-e 



Ci(x) • 



dt My* 

Eq. (22) is the equation of motion used in the simulation. 



( 20 ) 



( 21 ) 



( 22 ) 
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